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Abstract 

Many biological channels perform highly selective transport without direct input of metabolic energy and 
without transitions from a 'closed' to an 'open' state during transport. Mechanisms of selectivity of such 
channels serve as an inspiration for creation of artificial nano-molecular sorting devices and bio-sensors. To 
elucidate the transport mechanisms, it is important to understand the transport on the single molecule level 
in the experimentally relevant regime when multiple particles are crowded in the channel. In this paper we 
analyze the effects of inter-particle crowding on the non-equilibrium transport times through a finite-length 
channel by means of analytical theory and computer simulations. 

PACS numbers: 87.10.Ca, 87.10.Mn, 87.85.Rs 
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I. INTRODUCTION 



The functioning of living cells depends critically on molecular transport through various trans- 
port channels [1]. Many of them function without a direct input of metabolic energy and without 
a movable 'gate' that would involve transitions from an 'open' to a 'closed' state during transport. 
Nevertheless, such channels are selective, efficient and fast. Examples include porins, Nuclear 



Pore Complex and others [2,3,4, 5,161. The functioning of such channels has served as an inspira- 

:nnn^ 

tion for the creation of artificial biosensors and nano-molecular filters (7L 



a aim Hit iij, 



U that 



promise to play an ever increasing role in nano-technological and nano-medical applications, such 



:tion I71. 112 



and design of antibiotic drugs optimized for penetrating the cell nl5\\ . Such man-made channels 



110. 



also serve as testbeds for examining models of biological transport II 10 . 

Biological and artificial transport channels, such as those mentioned above, usually contain 
a passageway through which the molecules translocate by diffusion. From recent experimental 
and theoretical work, it has become increasingly clear that in many cases the transport selectivity 



of such channels is not dictated merely by molecule size, but is controllec 



by transient binding 



13, 



1811 . The crucial 



of the transported molecules inside the passageway 
insight into understanding the transport selectivity of such channels is that even in the absence 
of any physical barrier for the entrance to the channel, the probability of a particle to translocate 
through it is low (of an order of the aspect ratio of the channel) II 191]. T ransient trapping (due to 
binding) inside the channel overcomes this 'dimensionality barrier' [isl Jjl 2(], 21]. However, if 
the molecules are trapped in the channel for too long, the channel becomes crowded and transport 
is diminished. The interplay of these two effects provides a basis for selective transport, whereby 
only the molecules that are trapped in the channel for an optimal time transit through the channel 



with a high flux ||2 



0, 



21, 



22L 



230. Related mechanisms have been known in the context of carrier 

"(j, 24, 25]. Theoretical models that include 

,y,0, 



22, 



23 



3] provide 



assisted membrane transport as 'facilitated diffusion' yj 
the transient trapping combined with the effects of confinement |]16l . 
a good explanation of the behavior of the mean flux through nano-channels and show a good 
agreement with the experimental data lll7n . 

ogical perspective, transport of a single molecule can constitute a signifi- 
Thus, it is important to understand the transport through such channels 



However, from a bio 



cant signalling effect 

on the single molecule level. Advances in fluorescent microscopy and other methods allow one 



2 



to follow the transport of individual molecules through a channel H 121 1271. 128L |29l 130. 13 ill . Single 
molecule tracking experiments provide a wealth of information about the transport mechanisms, 
which is not accessible from the measurements of the bulk flux through the channel. The kinetics 
of transport of a single particle through the channel in the absence of other particles is well under- 



stood 119, 



32L 



33, 



3411 . In this paper, we analyze the effects of crowding of the particles inside the 



channel on the transport times of individual particles in the experimentally relevant regime when 
a non-equilibrium steady state flux passes through the channel. 



II. SINGLE PARTICLE 

Here, we briefly review the kinetics of a single particle passing through the channel in order 
to explain the methods employed herein. The channel is represented as a sequence of 'sites' 
1, ...,N. Inside, the particle performs diffusion-like random walk starting at the 'entrance' site 
1 and hopping between the internal sites 1 < i < N at an average rate r (for simplicity, we 
assume that the channel is uniform). The particle can leave the channel from the terminal sites 
1 and N with an average rate r Q . Transient trapping in the channel is described by choosing 
r Q < r. This hopping process is illustrated in Fig. 1. At any time t, the position of the particle 




r r r r 



FIG. 1: The channel is represented by a sequence of sites 1, N between which the particles can hop with 
rate r. The rate of hopping out of the channel from its ends is r D . In the single particle case, a particle starts 
at site 1 and hops inside the channel until it exits from either end. In the multi-particle case, the particles 
enter at site 1 with an average rate J, if its occupancy is less than the maximal allowed. The line shows the 
steady state concentration profile. 

in the channel is described by the vector of probabilities pi{t) to be at a particular site i: \p(t)) = 
(pi(t), ...pi(t)...p N (t)). We also define the vector \i) as a vector with the i'th element equal to 1 
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and all other elements equal to 0, so that (i\p(t)) = pi(t) (where (x\y) is the scalar product of 
the vectors \x) and \y)). The Master equation for the probability vector, describing the hopping 
through and out of the channel ends, can be written as (see the Appendix for details) 



j t \p(t)) = M.\p(t)). 
The formal solution of the equation (Q} can be written as \p(t)) 



the initial condition [33, 



(1) 



|p(0)), where \p(0)) is 



470; for a particle starting at site 1, |p(0)) = |1) = (1,0...0). The 



instantaneous probability flux to the right out of the channel is r pN(t), and the probability that 



the particle had exited the channel from the right side by time t is Pt> = J r p N 
The total probability to exit to the right P_> = P™ is 



220- 



r {N\e Mt \l)dt' = -r {N\M- x \l). 



(2) 



where (i\X\j) = Xij. After some algebra (see the Appendix for details), equation © gives for 
the total probability to exit from the channel on the right (the translocation probability): P^ = 
2+ (jv-i) r / r ' m accord with previous works QJ/ZL Ll9[|- Note that P^ increases as r„ diminishes. That 



is, trapping of the particle in the channel increases the translocation probability [|16 



19, 



3 



24, 



We now calculate the directional mean exit times. The probability distribution of the exit times 



to the right is 
to exit the right is 



1 d 

P_ dt 



1 - P* 



r oPN {t)/p^ m 



22, i 



t'f^{t')dt' = r {N\(M- 1 ) \1)/P^. 
Similarly, the mean first passage time to the left is 

^(0=r o (l|(M- 1 ) a |l>/P^. 



Thus, the mean time 



(3) 



The mean time to exit from any of the ends is 



(4) 



T 



* (Ptf(t) + Pi(t)) dt = T^P^ + T^P^. 



(5) 



Using the equations above, we obtain explicit expressions for the mean times: 

N (6P^ + P^ (N(N - 3) + 2) (r D /r) 2 ) 



1\ 



6r 

N (6P^P^ + Pi (N(2N - 3) + 1) (r /r) 2 ) 



Qr n P^ 



T 



N 
2^ 



(6) 
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in agreement with previous results obtained in the continuum limit H21L 13411 . Note that the mean 
trapping time T is linearly proportional to the channel length N. Surprisingly, the mean time for 
the particle to exit to the left TV also scales like iV for N 1, due to the possibility of large 
excursions into the channel before it returns to the left end. By contrast, the mean exit time to the 
right has two distinct regimes. For short channels, or strong trapping (Nr Q /r <C 1), TV ~ y~, 



while for long channels, or weak trapping, (Nr Q /r 3> 1), ~ ^- (see also Fig. 3). Physically, 
for strong trapping, the bottleneck for the exit to the right is the release from the channel end, 
while for long channels and weak trapping the exit time is dominated by the time it takes to 
diffuse through the channel from left to right. 

III. SINGLE PARTICLE ON THE BACKGROUND OF THE STEADY STATE FLUX 



When a finite flux J impinges onto the channel entrance, at any moment there can be many 
particles in the channel that might interfere with each other's passage and prevent the entrance 
of new ones. The particles in the channel obey the same kinetics as the single particles, with a 



condition that a site can contain up to a maximal number of particles m. Following II17 . 



37 



35, 3 



380, the system can be described in terms in terms of site occupancies = \n)^ s . For constant 



J, a non-equilibrium steady state is establishes 



be solved exactly: n\ 



JP^(l+(N-i)r /r) 
r +JP^/m 



iniM 



ie steady state profile of a uniform channel can 



36. 



37, 



4QD. 



We now turn to the main results of this paper - how does the crowding, when many particles are 
present in the channel, affect the transport times of individual particles within the non-equilibrium 
steady-state flux. To the best of our knowledge, no exact analytical solution exists in this case. 
The transport of an individual particle can be viewed as occurring on the background of the steady 
state density profile \n) ss . In the mean field approximation, the probability pAt) of a particle to be 



present at a given site is described by the following equations [38, 



39 



40, 



at m m 



(7) 



-rpi(l 



rpi(l 



m m 

with the appropriate boundary conditions (see the Appendix for details). Using matrix notations: 



j t \p{t)) = M-.\p{t)). 



(8) 



Explicit matrix elements of M ss are given in the Appendix. As in the single-particle case above, 
the linear equations (17181) can be solved analytically. 



To test the feasibility of the mean field approximation, we compared the probability of a particle 
to exit to the right, computed using the exact solution for the steady state density with the mean 
field result (see below). First, the average exit flux to the right is ,/_> = r n^, which yields for the 
probability of an individual particle within this steady state flux to exit to the right II 1711 : 

pss = J-* = 1 9 

J{l-nx/m) 2 + (N-l)r /r 

On the other hand, from the mean field approximation of eq. ([?]) Pff = — r (A|(M ss ) _1 |l) (see 
Eq. ©). Using the expressions for nf, after some algebra we get the same result as the exact 
expression, eq.©. Thus, the mean field approximation yields an exact result: the probability 
of an individual particle to exit to the right is not affected by crowding and is the same as in 
the single-particle case (at least for uniform channels) Q3VL |40J]. The directional mean exit times 
can be calculated by repeating the same algebra as for the case of a single particle, but with 
M ss instead of M. We find that the mean trapping time is T ss = - surprisingly, like the 
translocation probability, the mean trapping time is also not affected by the crowding. By contrast, 
the directional times to exit to the right and to the left, T*^, respectively, do change due to inter- 
particle interactions, compared to the single particle case. After some algebra (see the Appendix 
for details), one gets for the mean time to exit to the left: 



TZ =r (l\(M ss y 2 \l}/P< 
rm (JP^ + mr ) 2 



(10) 



J 3 r P^P^ 



rm 



JP- 



rm 



JP- 



Nm 

IppZ 



mr Q + JP^ + ^ 



where tp(x) 



dx 



r(x); T(x) is a 7-function. The mean exit time to the right can be obtained in a 



similar fashion, using equation ©. The dependence of the exit times on the impinging flux J is 
illustrated in Fig. 2. Unlike the exit probabilities, crowding increases the mean time to exit to the 
right T_>, and decreases the mean time to exit to the left TV. Interestingly, however, the qualitative 
dependence on the channel length N is similar to the single-particle case (J=0), as shown in Fig. 
3. Importantly, the transport times remain finite even in the fully jammed regime (J — > 00), when 
the flux through the channel saturates to its maximal value. In particular, the mean time to exit to 
the right tends to — — , while the mean time to exit to the left tends to — (see the Appendix for 
details). 

In order to corroborate the results of the mean field approximation and to investigate the limits 
of its validity, we performed computer simulations of the transport through the channel using 
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FIG. 2: Ratios of the mean exit times in the jammed regime to the single-particle times. Upper lines - 

T s _^/T^. Lower lines - T^/T< Solid red lines - analytical solution, dotted lines - simulations; N = 6, 

r = 1, r a = 0.1, m = 3 for all lines. Inset: Same for single file transport m = 1. 



a variant of the Kinetic Monte Carlo algorithm Il42h . Both the simulations and the mean field 
results show that the forward times are increased due to the jamming while the backward times 
are decreased. The increase in forward exit times is easily understood considering the reduction 
in hopping rates inside the channel due to crowding. The origin of the decrease in the backward 
exit time is more subtle: the crowding increases the number of particles which hop backwards out 
of the channel immediately after their entrance. For wide channels that can accommodate more 
than one particle at each site, the mean field results for the directional transport times and 
TV agree closely with the simulations. For strictly single file channels (m = 1), the mean field 
approximation underestimates the actual value of the exit time to the right TV and overestimates 
the exit time to the left TV, but still reproduces the right qualitative dependence of the times on 
the flux J and other parameters. The reason for the underestimation of the time to exit to the right 
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FIG. 3: Dependence of the transport time in the crowded regime on the channel length N. Dashed red lines 
- mean field results; symbols - corresponding simulations; r = 1, r = 0.01, J = 0.5. The dependence on 
N is qualitatively similar to the single particle case - see text. The inset shows the transition between small 
N (linear dependence) and large N (quadratic dependence) regimes; logarithmic scale; r = 1, r a = 0.01, 
J = 0.1. 



is that the mean field approximation neglects the correlation between successive jumps (a particle 
hopping to one of its neighbor sites leaves behind it a vacancy and thus has a higher probability 
to hop back to the same site in the next jump). Interestingly, the simulations show that the mean 
field result for the mean trapping time T is exact (at least for a uniform channel). The mean field 
approximation can be improved by taking into account the correlations in the jumping rates of 
the neighboring particles and the fluctuations of the density around its mean value. Correction 
to the mean field diffusion rate of a tracer particle in equilibrium conditions were calculated in 



[3 



43Q in the framework of effective medium theory. In general, such corrections to the mean 



field improve the approximation including our case (data not shown); however systematic analysis 
of such corrections lies outside the scope of the present work. 
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IV. SUMMARY 



To summarize, we have analyzed the effects of crowding and inter-particle competition for 
space on the transport times through narrow channels of finite length under a non-equilibrium 
steady state condition. The results of the mean-field analysis are corroborated by computer sim- 
ulations. We have shown that in uniform channels the jamming increases the forward exit time, 
while decreasing the backward exit time. Surprisingly, jamming does not affect the mean dwelling 
time in uniform channels. The situation might be different in non-uniform channels; however, the 
mean field approximation should provide a qualitatively correct picture even in this case [38], full 
discussion of which lies beyond the scope of the present work and will be discussed elsewhere. 
The model provides a theoretical framework for analysis of single molecule transport through bio- 
logical and artificial nano-channels. The parameters of the model, the rates r Q and r, can be related 
to the experimentally controlled factors such as diffusion coefficients inside and outside the chan- 
nel and the binding affinity of the molecule in the channel. It is also important to emphasize the 
difference between the results of this paper and the well studied case of tracer diffusion in infinite 

n 

single-file channels [44]. Finally, we note that the methods of this work can be extended to treat 



460, confor- 



arbitrary molecular signalling pathways, such as multi-step enzymatic reactions H45L 
mational transitions of ion channels [47] and other systems |48J]. 

The authors are thankful to R. Groger, I. Nemenman, B. Munsky, A. Perelson, K. Rasmussen, N. 
Sinitsyn for stimulating discussions. This research was performed under the auspices of the U.S. 
Department of Energy. 



APPENDIX: DETAILES OF THE CALCULATIONS IN THE MAIN TEXT 



1. Single particle 



At any time t, the position of the particle in the channel is described by the vector of probabil- 
ities pi(t) to be at a particular site i: \p(t)} = (p-i (t )... .Pi (t) ...p N (t) ) , so that {i\p(t)) = Pi(t) The 



probabilities pi(t) obey the following equations [|2lL 



32, 



33] 



^-Pi(t) = r(pi_i + p l+1 - 2pi) for 1 < % < N 
at 



(A.l) 



with the boundary conditions 

d 



dtP 1 ^) = ~( r ° + r )P! + r ^ 2 ) and di PN ^ = ~( r ° + r )P N + r PN-i- (A.2) 



Equations (IA.1IA.2I) can be written in a matrix form as 

j t \p(t))=M-\p(t)), (A3) 

with 

M ifi = -2r and M i>i±1 = r for 1 < i < N, (A.4) 

and 

M hl = -r-r Q - M NtN = -r-r ; M lj2 = r; M^iv-i = r. (A.5) 
a. Explicit solution of single particle equations in terms of matrix elements 

Here we re-derive the solutions obtained in the main text, using the standard methods of linear 



algebra [49]. Assume an arbitrary Markov process that can be in N states (such as defined by 
equation (3)). Time evolution of its probability distribution p(t) = (pi,P2, —,Pn) an d can be 
described by the following matrix equation (the equation (3) is an example): 

p = U -p (A.6) 

where p is an A^-dimensional vector, and U is an N x N matrix. Let us denote the eigenvalues 
of the matrix U as uj 1 ...uj i ...uj N and the corresponding eigenvectors as v 1 , ...v N . Then the 
general solution is 

N 

p{t) = ajiPePt*, (A.7) 



3=1 

where ai.-.a^ is a set of numerical coefficients. In other words, 



N 

ft (t) = ^aX'e^. (A.8) 

3=1 

The coefficients aj can be determined from the initial condition: 

TV 

Pl (0) = J2ajvl (A.9) 

3=1 

which can be written as 

p(Q) = V.a J (A. 10) 
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where 



v 



v' 



so that Vij = v\ . Finally 



and 



a = y- x ■ p(0) 



(A.11) 



JV 



(A.12) 



For the initial condition Pi(0) = (^i, we get ak = (V -1 )^. Now, the matric U is diagonalized 
to its diagonal form 

/ w, ... 0\ 



w 



... Wi ... 



VO ... w N ) 



by the transformation W = V^UV, or equivalently U = VWV' 1 04911 . 
Thus, the probability to be in state i at time t is 



N 



N 



Pi (t) = 5>X'exp^ = ^V" 1 ^^ = (^V- 1 )^ = (e u \ (A.13) 
i=i i=i 

and in particular, 

P*(*) = {e Ut ) m (A.14) 
The probability flow to exit to the right is r Q p N (t), and the total probability of exit to the right 

is 



P- = r, 



r°° N i 

/ p N (t)dt = -r Q a j v N — 

* i 



(A.15) 
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in agreement with equation (2) in the main text. This result can be also obtained using the 
following reasoning. Instead of considering a single particle hopping through the states, starting 
at the position 1, let us consider the steady state where a flux J enters to a position 1, with a 
steady state probability distribution p. Then the probability to exit to the right is the ratio of the 
transmitted flux to the entrance flux: r Q p N / J. We have 

= U-$+J (A.16) 

and therefore p = —U~ l ■ J so that p N = —Uz\ because Jj = JSn. 



The probability distribution of exit times to the right is simply r pjsi{t) (12111 and any moment of 
it can be calculated easily. For instance, the mean first passage time to exit to the right is: 

/•oo n 

f^ = r tP N {t)dt = ro^iV-^jiVNj— = r (V^^r^ = r (U 2 ) N \ (A.17) 
Jo j=1 ™j 

in agreement with the equation (3) in the main text. 



2. Steady state 

In the case of current J impinging on the channel entrance, one can describe the system in 
terms of average site occupancies n,;, whose kinetics is described in the mean field approximation 



by the following equations II17 . 



35 



31 



^-rii = rn i _ 1 (l - — ) + rra i+1 (l - —) 
at mm 



/-, n>i-l\ n i+l 

rndl ) — rrii{l 



m 



r(rij_! + n i+ i — 2rij) for 1 < i < N. 



(A. 18) 



where m is the maximal site occupancy. The boundary conditions at sites 1 and N are 



d_ 

dt 
d 

Tt 



dt ni 



n N 



— Or + r )n\ + rn 2 + J(l ] 

m 

-Or + r )n N + rn^-i- 



(A.19) 



In a matrix form: 



d_ 

dt 



n(t)) = M J ■ \n(t)) + J 



(A.20) 



where the matrix M J is the same as M with the only change M( x = —J/m — r — r and J 



(J, 0, ...0). Note that for an internally uniform c 



field equations (Eqs. (IA.18IA.20I) ) are exact 1135 . 



lannel (as the one described in Fig. 1) the mean- 



37D 



12 



The steady state density profile can be obtained from Eq. (1A.20I) as \n) ss = — (^M J ^j ■ J, or 
more specifically as: 

J(l+ (N-i)*-) 

1 r (2 + (N-iy f ) + i(i + (N-iy f y 1 ■ } 

The average exit flux to the right is J^, = r n s ^. This together with Eq. (IA.21I) yield the probability 
of an individual particle within the flux to exit to the right: 

p ss — _ = \ (A 22) 

J{l-n x /m) 2 + (N- l)r a /r' 

As already established before, the exit probability of individual particles to exit to the right is the 

same as in the single-particle case (at least for uniform channels), even though they are interfering 

with each other's passage through the channel lv\ . 

However, crowding does influence transport and is manifested in obstruction of the entrance 

site. The transport efficiency, defined as the ratio of the exit flux to the right to the total 

impinging flux J, Eff^ = ^f-, decreases with J due to jamming at the entrance. 

Eff — = r _2_ (A 23) 

J 2r + (N - l)r Q + J(l + (N - l)r /r)/m 

3. Derivation of the analytical expressions for the mean exit times in the jammed regime 

The boundary conditions of equation (8) of the main text, describing the probability of the 
tagged particle are: 

-f-pi = — r Q pi — rpi(l — nf/m) + rp 2 (l — n\ s Im) 
at 

-^Pn = -r a PN - rp N (l - n s x_ x /m) + rp N -i{l -n^/m) (A.24) 

Using the matrix form of Eq. (9) of the paper, the elements of the matrix M ss are given by 

M?? = -r(2 - ntJm - n^Jm) = -2r(l - nf/m); (A.25) 
and M- - ±1 = r(l - nf/m) for 1 < i < N, 



and 



M{\ = -r(l - -4 - r ; M S J N = -r(l - - r Q ; (A.26) 



m ' m 

ry\ SS ry-, SS 

M% =r(l-^); M- i7V _i=-(l- — 

m m 
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The average (over particles actually exited to the left) time to exit to the left is 

TZ=r (l\{(M s T 1 ) 2 \l)/P^ 

where 

l + (iV-l)r /r 
~* 2 + (JV-l)r /r" 

In order to obtain an explicit expression for we define 

\W) = D{M SS )' 1 |1) 

and 

(Q\ =D(l\(M ss )-\ 

In the equations above we introduced the notation for the determinant of M ss , D = det 
The elements of these vectors are given by 

N-l 

W n = r N ~ 2 {A N r + nB N r + {N - n) r ) J] (A N + kB N ) 

k=2 

and 

N-l 

11 (A N + kB N ) 

,N-2 (A „ i „ D „ i / AT „\ „ \ fc=l 



Q n = r JV " 2 (A^r + nB N r + (JV - n) r Q ) 



A N + nB N 
In the above expressions 

r Q (2r + (JV - 1) r D - J/m) 



A 



N 



(J/m) (r + (JV - 1) r ) + r Q (2r + (JV - 1) r ) 
and 

Jr /m 

-DAT 



(J/m) (r + (JV - 1) r ) + r Q (2r + (JV - 1) r ) ' 
The mean escape time to the left is then 

N 



n=l 



Using the notations above we can express D as: 

N-l 

D = r N ~ 2 r (2A N r + (JV + 1) B N r + (JV - 1) r ) ] [ (A N + kB N ) . 



k=2 
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Substituting D into the expression for T ^ one obtains 

JV 



(An + B N )J2 iAiE±I ^- 



r+{N-n)r f 



-nB A 

TZP^ = — — 7TT~"T^ : — ttt — 7— ^i- ( A - 37 ) 



r (2A N r + (JV + 1) + (JV - 1) r ) 
Performing the summation, we get for the average time to exit the channel to the left 

TZP^ (N(B N r - r )(2A N (B N r + r ) + B N (B N (N + l)r + (3JV - l)r ))) 



(Ajv + B N ) 2B 2 N r Q {2A N r + (JV + 1) i^r + (iV - 1) r 

r^Ajv + JVJ^) 2 (V f|f + N + l) - t/f ( 

+ V 7/ . (A.38) 

B% (2A N r + (JV + 1) B N r + (JV - 1) r f 



A N +B N 
B N 



Here ^ (x) = dln ^ z ^ | z=x , where is the 7-function. 

Substituting the expressions for A N and B N , we get the explicit expression for the average time 

as 

- ss _ JV (2r p ( 2r + r (iV - 1)) + £ (4r + 3r D (iV - 1))) 



2 (J/m) 2 (r + r (JV - 1)) 

(J/m) (r + r (JV - 1)) + r (2r + r Q (JV - 1)) / / at 2r + (JV - 1) r \ ( / 2r + (JV — 1) 

ip JV H — J — ^ 



( J/m) 3 r a (r + r G (JV — 1))) \ \ J/m J \ J/m 

(A.39) 

In the main text we use the definitions of the probabilities to exit to the right/left in order to simplify 
this cumbersome expression. In the single particle limit, J — > 0, the expression for TZ reduces 
to the previously obtained single particle expression (Eq. (6) in the paper). At the other extreme 
when the input flux J — > 00, the mean escape time is 

lim TZ = ; + (i "- 1)r ; , = ^. (A.40) 
J-00 - r (2r + (JV - l)r ) r 

The mean time to exit to the right and the mean trapping time can be obtained in a similar fashion. 
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